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Abstract 

This paper considers the quantization problem on the Grassmann manifold Q n , P , the set of all p-dimensional 
planes (through the origin) in the n-dimensional Euclidean space. The chief result is a closed-form formula for 
the volume of a metric ball in the Grassmann manifold when the radius is sufficiently small. This volume formula 
holds for Grassmann manifolds with arbitrary dimension n and p, while previous results pertained only to p = 1, 
or a fixed p with asymptotically large n. Based on this result, several quantization bounds are derived for sphere 
packing and rate distortion tradeoff. We establish asymptotically equivalent lower and upper bounds for the rate 
distortion tradeoff. Since the upper bound is derived by constructing random codes, this result implies that the 
random codes are asymptotically optimal. The above results are also extended to the more general case, in which 
G n ,q is quantized through a code in G n ,p> where p and q are not necessarily the same. Finally, we discuss some 
applications of the derived results to multi-antenna communication systems. 

Index Terms 

the Grassmann manifold, rate distortion tradeoff, MIMO communications 

I. Introduction 

The Grassmann manifold Q n:P (L) is the set of all p-dimensional planes (through the origin) in the 
n-dimensional Euclidean space L n , where L is either R or C. It forms a compact Riemann manifold of 
real dimension (5p (n — p), where (5 = 1 when L = R and (5 = 2 when L = C. The Grassmann manifold 
provides a useful analysis tool for multi-antenna communications (also known as multiple-input multiple- 
output (MIMO) communication systems). For non-coherent MIMO systems, sphere packings of Q n ^ (L) 
can be viewed as a generalization of spherical codes [l]-[3]. For MIMO systems with partial channel 
state information at the transmitter (CSIT), which is obtained by finite-rate channel-state feedback, the 
quantization of beamforming matrices is related to the quantization on the Grassmann manifold [4]-[6]. 

The basic quantization problems addressed in this paper are the sphere packing bounds and rate distortion 
tradeoff. Roughly speaking, a quantization is a representation of a source in Q np (L). In particular, it maps 
an element in Q n p (L) into a subset of Q np (L), known as a code C. Define the minimum distance of a 
code 5 = 5 (C) as the minimum distance between any two codewords in the code C. A sphere packing 
bound relates the size of a code and a given minimum distance S. Rate distortion tradeoff is another 
important aspect of the quantization problem. A distortion metric is a mapping from the set of element 
pairs in Q np (L) into the set of non-negative real numbers. Given a source distribution and a distortion 
metric, the rate distortion tradeoff is described by the minimum expected distortion achievable for a given 
code size, or equivalently the minimum code size required to achieve a particular expected distortion. 

There are several papers addressing the quantization problem for Grassmann manifolds. In [7], an 
isometric embedding of Q niP (R) into a sphere in Euclidean space ]R2( m_1 )( m+2 ) is given. Then, using the 
Rankin bound in Euclidean space, the Rankin bound in Q n ^ p (R) is obtained. Unfortunately, this bound 
is not tight when the code size is large. Instead of resorting to an isometric embedding, sphere packing 
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bounds can also be derived from analysis in the Grassmann manifold directly. Let B (<5) denote a metric 
ball of radius 5 in Q np (L). The sphere packing bounds can be derived from the volume of B (5) [3]. The 
exact volume formula for a B (5) in Q np (L) with p — 1 and L = C is derived in [4]. An asymptotic volume 
formula for a B (5) in Q ntP (L), where p > 1 is fixed and n approaches infinity, is derived in [3]. Based 
on those volume formulas, the corresponding sphere packing bounds are developed in [3], [5]. Besides 
the sphere packing bounds, the rate distortion tradeoff is also treated in [8], where approximations to the 
distortion rate function are derived via the sphere packing bounds on the Grassmann manifold. However, 
the derived approximations are based on the volume formulas in [3], [4] which are only valid for some 
special choices of n and p: either p = 1 or fixed p > 1 with asymptotic large n. 

The main contribution of this paper is to derive a closed-form formula for the volume of a small ball 
in the Grassmann manifold. Based on this formula, sphere packing bounds are derived and rate distortion 
tradeoff are accurately quantified. Specifically: 

1) An explicit volume formula for a metric ball in Q np (L) is derived when the radius 5 is sufficiently 
small. It holds for Grassmann manifolds with arbitrary dimensions while previous results are only 
valid for either p = 1 or a fixed p with asymptotically large n. The main order term of the volume 
is c n) p )P) /3(5^ p< - n ~ p ' ) for a constant c n:P:Pj p depending on n, p and (3. Lower and upper bounds on the 
volume formula are also derived. 

2) Based on the volume formula, the Gilbert- Varshamov and Hamming bounds for sphere packings are 
obtained. For the distortion rate function, a lower bound is established via sphere packing argument 
and an upper bound is derived via random-code argument. The bounds are in fact asymptotically 
identical, and so precisely quantify the asymptotic rate distortion tradeoff. Since the upper bound 
is actually derived from the average distortion of random codes, it follows that random codes are 
asymptotically optimal. 

3) The volume formula and the results on the rate distortion tradeoff are extended to a more general 
plane matching problem. In this plane matching problem, a plane from the code C C Q n;p (L) is 
chosen to match a random plane Q £ Q nA (L) to minimize the distortion, where p and q are not 
necessarily the same. For plane matching, a metric ball in Q n q (L) centered at a plane in Q np (L) 
is studied. The volume formula is derived for such a ball with sufficiently small radius. The rate 
distortion tradeoff is also quantified by the same method as above. 

4) As an application of the derived quantization bounds, the information rate of a MIMO system with 
finite-rate channel-state feedback and power on/off strategy is accurately quantified for the first time. 
Since the corresponding Grassmann manifold for most practical MIMO systems has p > 1 and small 
n, the quantization bounds derived in this paper are necessary. 

The paper is organized as follows. Section |TT| provides some preliminaries on the Grassmann manifold. 
Section HTD derives the explicit volume formula for a metric ball in the Grassmann manifold. The corre- 
sponding sphere packing bounds are obtained and the rate distortion tradeoff is accurately quantified in 
Section |TVJ An application of the quantization bounds to MIMO systems with finite-rate channel-state 
feedback is detailed in Section |V| Section |VT| contains the conclusions. 

II. Preliminaries 

This section presents a brief introduction to the Grassmann manifold. A metric and a measure on the 
Grassmann manifold are defined, and the problems relevant to quantization on the Grassmann manifold are 
formulated. For completeness, we also extend the quantization problem to a more general plane matching 
problem. 

A. Metric and Measure on Q n ^ p (L) 

For the sake of applications [4]-[6], the projection Frobenius metric {chorda! distance) is employed 
throughout the paper although the corresponding analysis is also applicable to the geodesic metric [3]. 
For any two planes P,Q£ Q n , P QL), we define the principle angles and the chordal distance between P 
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and Q as follows. Let u x G P and V] G Q be the unit vectors such that 



ujvi 



is maximal. Inductively, 



let Uj G P and Vj G Q be the unit vectors such that u]uj = and vjv,- = for all 1 < j < i and 



is maximal. The principle angles are defined as Q { = arccos 
distance between P and Q is given by 



for i = 1, • • • ,p [7], [9]. The chordal 



d c (P,Q) 



]Tsm 2 ^. (1) 
^ i=i 

The invariant measure /i on Q n ^ (L) is defined as follows. Let O (n) and U (n) be the groups of n x n 
orthogonal and unitary matrices respectively. Let A,B 6 O (n) when L = R, or A, B G £/ (n) when 
L = C. For any measurable set M C ^„ iP (L) and arbitrary A and B, 

(AM) = n (M) =fi(MB). 

The invariant measure defines the uniform/isotropic distribution on Q n , p (L) as well [9]. 

B. Quantization on Q n<p (L) 

Given both a metric and a measure on Q n>p (L), a quantization on the Grassmann manifold can be well 
defined. Let C be a finite size discrete subset of ^„ iP (L). A quantization is a mapping from the Q njP (L) 
to the set C (also known as a code), i.e., 

q:£ n , p (L)->C. 

An element in the code C is called a codeword. Thus, roughly speaking, a quantization is to use a subset 
°f Qn,p (L) to represent the whole space. 

Sphere packing bounds relate the size of the code to the minimum distance among the codewords. Let 
5 be the minimum distance between any two codewords of a code C and B (8) be a metric ball of radius 
8 in the Q TliP (L). If is any positive integer such that Kfi (B (5)) < 1, then there exists a code C of size 
K + 1 with minimum distance 5. This principle is called as the Gilbert-Varshamov lower bound, 

H{B(8)) 

On the other hand, |C| /x (P (5/2)) < 1 for any code C. The Hamming upper bound captures this fact as 

For more information about the sphere packing bounds, see [3]. 

Rate distortion tradeoff is another important aspect of the quantization problem. A distortion metric is 
a mapping, 

: £„, P (L)xC^ [0, +oo) , 

from the set of the element pairs in Q n p (L) and C into the set of non-negative real numbers. Throughout 
this paper, we define the distortion metric as the square of the chordal distance, d (■, ■) = d\ (•, ■). Assume 
that a source Q is randomly distributed in Q 1hP {h). The distortion associated with a quantization q is 
defined as 

P^E[c) (Q,q(Q))]=E[d 2 c (Q,q(Q))}. 

The rate distortion tradeoff can be described by the infimum achievable distortion given a code size, 
which is called the distortion rate function, or equivalently the infimum code size required to achieve a 
particular distortion, which is called the rate distortion function. In this paper, the source Q is assumed to 
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be uniformly distributed in Q n>p (L). For a given code C C G n ,p (L), the optimal quantization to minimize 
the distortion is given by 1 

q (Q) = arg mm d c (P, Q) . 



The distortion associated with this quantization is 

D{C) = E 



min d\ (P, Q) 

Pec 



For a given code size K where K is a positive integer, the distortion rate function is then given by 2 



D* (K) = inf D (C) . 

C:\C\=K 



The rate distortion function is given by 



K* (D) = inf ICI 

D(C)<D 



(4) 



(5) 



C. An Extension: Plane Matching Problem 

For the sake of completeness, we extend the quantization problem to a more general plane matching 
problem. The plane matching problem involves planes from different spaces Q ntP (L) and Q n>q (L) where 
p and q are not necessarily the same. 

To formulate the plane matching problem, we need to define the chordal distance between P G Q n>p (L) 
and Q 6 Q nA (L). Without loss of generality, we assume that p < q. Using the same procedure described 



< P < 



Based on the 



in Section Pl-AL we are able to define the principle angles < 6>x < „ ,,. 

principle angles, the chordal distance between P E Gn, P (L) and Q G Q n , q (L) are defined as d c (P, Q) = 



Y^=i sm2 ®i- I n m i s wa y 5 the definition of chordal distance in ([TJ is just a particular case of the general 
definition. 

Now consider the plane matching problem. Intuitively, the plane matching problem is to choose a plane 
from the code C C (L) to match a random plane Q G Q n , q (L) such that the average distortion is 
minimized, where 1 < p < n and 1 < q < n are not necessarily the same. Formally, a plane matching is 
a map from the whole space of Grassmann manifold, e.g., Q nA (L), to the code C C Q n ^ v (L), 



q:g n>q (L)-C, 



such that 



D (C) 4 Eq 



min dl (P, Q) 

Pec cv 



is minimized. According to the same principles in the quantization problem, the rate distortion tradeoff 
can be extended to the plane matching problem. 

The ties, i.e. the case that 3-Pi, Pi £ C such that d c (Pi, Q) = min d c (P, Q) = d c (P2, Q), are broken arbitrarily as they occur with 
probability zero. 

2 The standard definition of the distortion rate function involves the code rate, which is log 2 K. The definition in this paper is equivalent 
to the standard one. 
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III. Metric Balls in the Grassmann Manifold 

In this section, an explicit volume formula for a metric ball B (5) in the Grassmann manifold is derived. 
It is the essential tool to quantify the rate distortion tradeoff in Section [TV] 

The volume calculation depends on the relationship between the measure and the metric defined on 
the Grassmann manifold. This paper focuses on the invariant measure /i, which corresponds to the 
uniform/isotropic distribution, and the chordal distance d c . For any given P E Q n ^ v (L) and Q E Q n>q (L), 
define 

B P (5) = [QE Q n>q (L) : d c (P, Qj < 5 

and 

B Q (6) = IP E g rhP (L) : d c (P, q)<s). 



For the invariant measure fi, it has been shown that /x (Bp (5)) — fj, (Bq (5)) and the value is independent 
of the choice of the center [9]. It is convenient to denote B P (5) and Bq (5) by B (5) without distinguishing 
them. Then, the volume of a metric ball B (5) is given by 



H(B(5)) 



d/j e 



(6) 



where l<#i<f,---,l<0j,<f are the principle angles and the differential form d/ig is the joint 
density of the 6>j's, which is given in [9]— [1 1] and as well in Appendix IXI below. 

The following theorem calculates the volume formula and expresses it as an exponentiation of the 
radius. 

Theorem 1: When 8 < 1, the volume of a metric ball B (5) is given by 



PL [B (5)) = ( 1 + c« ^ + o (5 2 ) 



(7) 



where 



1 if L = R 

2 ifL = C 



c 



n,p,q,P 



r( 


§p(n- 










r( 


§p(n- 






1 




r( 


fp(n- 






1 




r( 


f p(n- 





and 



^n,p,q,f3 



„ r(|(n-i+l)) 

w=i r(|( 9 -i+i)) 

TTP r(f(n-i+l)) 

TT 1 Ai=i r (f(n-p-i+i)) 

n _ g r(f(n-i+l)) 

i=1 r(f(g-i+l)) 

n _ ? r(|(n-i+l)) 
i=1 r(f(n-p-i+l)) 



if 1 < p < q < f 

ifl<p<|<g<-n, and p + q < n 
iil<p<^<q<n and p + q > n 
ii^<p<q<n 



(8) 



1 



f P (n - g) 
|p (n - q) + 1 



(9) 



Proof: See Appendix IAI 
The following corollary gives the two cases where the volume formula becomes exact. 
Corollary 1: When 5 < 1, in either of the following two cases, 

1) L = C and q = p; 

2) L = R and q = p+ 1, 

the volume of a metric ball B (5) can be exactly calculated by 

n(B(8)) = c niPiq ,p5M n -«\ 

where c nyPyq ^ is defined in (HJ). 
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We also have the general bounds: 

Corollary 2: Assume 8 < 1. If L = R and p = q , the volume of B (8) is bounded by 



c n , P , P ,i8 p(n - p) <fi(B (8)) < c ml 8 p ^ (1 - 8 2 ) 



For all other cases, 



(1 - 5 2 ) |p(9 " P+1) - p c n , p ,^^(-^ <p(B (8)) < cw^-d. 
Proof: Corollary [T] and |3 follow the proof of Theorem [l] by tracking the higher order terms. ■ 
Theorem [T] is of course consistent with the previous results in [4] and [3], which pertain to special 
choices of n and p and are stated below as examples. 

Example 1: Consider the volume formula for a B (8) where p = q = 1. Without normalization, the total 
volume of Q^ x (C) is 2-n n / (n - 1)! and the volume of the B (8) is 2vr n 5 2 ( n - 1 V (n - 1)! [4]. Therefore, 

agreeing with Theorem [T] where (3 = 2 and c n 1 1 2 = 1. 

Example 2: For the case that p = q are fixed and n — > +00, an asymptotic volume formula for a B (5) 
is derived by Barg [3], which reads 

/ X \ /3«P+o(n) 

l*{B(S))= — . (10) 



VpJ 

On the other hand, asymptotic analysis from Theorem [T] gives 

{3p(n—q)+0 (log n) 

H(B(6)) 



/ x \ PP(n-q)+u(Logri) f fa \ \ 

(±) (i-(f(,-P + i)-i)^ + (^)) 



for asymptotically large n and any fixed 1 < p < q. The derivation follows the Stirling's approximation 
applied to — ^^^y log (c^p ^). In this setting, Theorem[l]is consistent with (ITOb and provides refinement. 

Importantly though, Theorem [T] is distinct from the above results in that it holds for arbitrary p, q and 
n. For a metric ball with parameter n not asymptotically large, i.e., p and q are comparable to n, it is not 
appropriate to use (TTOb to estimate the volume. A trivial example is that the p = q = n case. If p = q = n, 
the exact volume of B (8) for V5 > is the constant 1. The formula in Theorem [l] gives c n>n>n ^ = 1 

and Cn,n,n,p&^ n ~ q ^ — 1- However, the approximation (8 / \fnf n (formula (fTOb ) will give a small number 
much less than 1 when 8 is small. 

For engineering purposes, it may be satisfactory to approximate the volume of a metric ball B (8) by 
c n ,p,g,/38 l3p ^ n ~ 9 ' when the radius 8 is relatively small. Fig. [T] compares the exact volume of a metric ball © 
and the volume approximation c n>p>qj p8^ n ~ q \ In the simulations, we always assume p = q. The volume 
approximation becomes exact for the complex Grassmann manifold when 8 < 1. To calculate the exact 
volume without appealing to Theorem [H Monte Carlo simulation is employed to evaluate the complicated 
integrals in ©. Since 

fi(B(8)) = Pr{Q: d c (P,Q)<8} 

where P E Q ntP (L) is chosen arbitrarily and Q is uniformly distributed in the Q n:P (L), simulating the 
event {Q : d c (P,Q) < 8} gives fx(B(8)). The simulation results for the real and complex Grassmann 
manifolds are presented in Fig.[TJa) and[TJb) respectively. Simulations show that the volume approximation 
(solid lines) is close to the exact volume (circles) when the radius of the metric ball is not large. We also 
compare our approximation with Barg's approximation (8/^/pj np (from (fTOl ) for the case n = 10 and 
p = 2. Simulations show that the exact volume and Barg's approximation (dash-dot lines) may not be of 
the same order while the approximation in this paper is much more accurate. 
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IV. Quantization Bounds 

This section derives the sphere packing bounds and quantifies the rate distortion tradeoff for both 
quantization problem and plane matching problem. The results developed hold for Grassmann manifolds 
with arbitrary n, p and q. 

A. Sphere Packing Bounds 

The Gilbert- Varshamov and Hamming bounds for Q n v (L) are given in the following corollary. 
Corollary 3: When 5 is sufficiently small (5 < 1 necessarily), there exists a code C in the Q n ^ v (L) with 
size K and the minimum distance 5 such that 

(1 + 0(5*)) <K. 

For any code with the minimum distance S, 

*<C(j) (i + oO* 2 ))- 

Proof: The corollary is proved by substituting the volume formula Q into © and ©. ■ 
Remark 1: Applying Corollary El would provide sharper information on the higher order term. But we 
omit this. 

B. Quantization: the Rate Distortion Tradeoff 

The rate distortion tradeoff for quantization is characterized in this subsection 3 . Here, we assume that 
the quantization is on Q ntP (L), the source is uniformly distributed in Q UjP (L) and the distortion metric 
is defined as the square of the chordal distance. The derivation is based on the volume formula in Q. A 
lower bound and an upper bound on the distortion rate function are established. Denote the size of the 
code by K. Then the lower and upper bounds are asymptotically identical when p is fixed, n and the 
code rate log 2 K approach to infinity with a fixed ratio. Therefore, these bounds precisely quantify the 
asymptotic rate distortion tradeoff. Note that the upper bound is the average distortion of random codes. 
Random codes are asymptotically optimal. 

The following theorem gives a lower bound and an upper bound on the distortion rate function. 

Theorem 2: When K is sufficiently large {{c n ^^ p ^K)~ ' m™-p) < 1 necessarily), the distortion rate 
function is bounded as in 

Pp (n — p) 



(3p (n — p) + 2 

2r 



(c n)PiP ^Krm^ (l + o (1)) < D* (K) 



Remark 2: For engineering purposes, the main order terms in dTTb are usually accurate enough to 
characterize the distortion rate function. The details of the (1 + o (1)) correction are spelled out in Theorem 
|3] below, where a quantization is viewed as a specific case of plane matching. 

The lower bound and the upper bound are proved in Appendix |H] and O respectively. We sketch the 
proof as follows. 

The lower bound is proved by a sphere packing argument. The key is to construct an ideal quantizer, 
which may not exist, to minimize the distortion. Suppose that there exists K metric balls of the same 
radius 5o packing and covering the whole Q n ^ p (L) at the same time. Then the quantizer which maps each 
of those balls into its center gives the minimum distortion among all quantizers. Of course such an ideal 
packing may not exist. It provides a lower bound on the distortion rate function. 



3 For compositional clarity, the results for plane matching is summarized in a separate subsection IIV-CI 
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The basic idea behind the upper bound is that the distortion of any particular code is an upper bound 
of the distortion rate function and so is the average distortion of an ensemble of codes. Toward the 
proof, the ensemble of random codes C ran d = {Pi, • ■ • , Pk} are employed, where the codewords P;'s are 
independently drawn from the uniform distribution on Q UjP (L). For any given code C ran d, the corresponding 
distortion is given by 

D (Cand) = Eq 



min d c (Pi, Q) 

P PC j 

where Q 6 Q n>p (L) is a uniformly distributed plane. It is clear that D* (K) < Ec rand 
to calculate E Crand [D (C rand )]. Note that 



[D (C rand )]. We want 



[D(C, 



rand J 



E, 



min d\ (P, Q) 

P PC j 

± 2 errand 

min d 2 c (P u Q) 

% tLrand 



[min d 2 c (Pi, Q)\ should be inde- 



Since C ran d is randomly generated from the uniform distribution, Ec r . 
pendent of the choice of Q. Therefore, 

E Crand [D (Cr an d)] = E Crand [min d\ (P, Q)] 

for any fixed Q. By the volume formula and the extreme order statistics, we are able to calculate the 
distribution of d\ (Pi,Q) (for Vi) and min d\ (Pi,Q). In appendix O we prove that for any given Q 6 
Q n>p {h), K~i ■ E Crand [D (C ra nd)] converges to a constant as K approaches infinity. Therefore, an upper 
bound of the distortion rate function is obtained for asymptotically large K. 

The rate distortion function is directly related to the distortion rate function. The following corollary 
quantifies the rate distortion function. 

Corollary 4: When the required distortion D is sufficiently small (D < 1 necessarily), the rate distortion 
function satisfies the following bounds, 



f3p(n-p) 



Pp (n — p) 



Cn,p, P ,/3 \ 2T 



< 



D 



;i + (i))<ir(p>) 



(3p(n-p) 

1 ff3p(n-p)+2 



D 



0p(n—p) 
2 



(1 + 0(1)). 



(12) 



c n ,p,p,p \ (3p(n-p) 

To investigate the difference between the lower and upper bounds in (frill , proceed as follows. Since 
the exponential terms are the same in both bounds, focus on the coefficients. The difference between the 
two bounds depends on the number of real dimensions (3p (n — p) of the underlying Grassmann manifold. 
There are three cases to consider. 

Case 1: (3p (n — p) = 0. This only occurs n = p. Then the whole Q n ^ n (L) contains only one element 
and no quantization is needed essentially. 

Case 2: (3p(n — p) = 1. This happens if and only if L = R, n = 2 and p = 1. In this case, it can be 
verified that the principle angle 9 between a uniformly distributed Q G £2,1 (K.) and any fixed P e Q 2 ,i (^) 
is uniformly distributed in [0, |] . From here, the optimal quantization can be explicitly constructed. Since 
there exists K metric balls with radius sin such that those balls not only pack but also cover the whole 
£2,1 (R.)> the quantizer mapping those balls into its center is optimal. The distortion rate function can be 
explicitly calculated as 

Case 3: f3p (n — p) > 2. For this general case, an elementary calculation shows that 

(3p(n-p) 2 _ ( 2 



1 

2 ~ f3p (n - p) + 2 



< r 

(3p (n — p) \ (3p (n — p) 



< 1, 
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and we expect the difference between the two bounds to decrease as n approaches infinity. Indeed, the 
following corollary shows that the lower and upper bounds are asymptotically the same. 

Corollary 5: Suppose that p is fixed, n and the code rate log 2 K approach to infinity simultaneously 
with f = lim log2 K . If the normalized code rate f is sufficiently large (p2~~Pp r < 1 necessarily), then 

(n,K)— »+oo 

lim D* (K) = p2~& . 

On the other hand, if the required distortion D is sufficiently small (D < 1 necessarily), then the minimum 
code size required to achieve that distortion satisfies 

lim log 2 ^P) = ^ lQg2 (P_). 
{ritK )^+oo n 2 2 \DJ 

Proof: The leading order is read off from 

PP (n ~ V) , r - 

lim — ; : = 1 = lim — ; -1 



n^+oo f3p (n — p) + 2 rw+oo f3p (n — p) \(3p(n — p) 

2 

lim {cn,p, p ,p) mn ~ p) = P, 



n— »+oo 



2 



and 



lim K v(n--v) =2 Pp 



K-C n # tPtj3 — > oo. 



That the (1 + o (1)) multiplicative errors fall into place is the content of Theorem ■ 
The lower and upper bounds asymptotically agreeing accurately quantifies the distortion rate function. 

Since the upper bound is actually derived from the average distortion of random codes, this implies that 

random codes are asymptotically optimal. 

As a comparison, we cite the distortion rate function approximation derived in [8]. For Q n , p (C) with 

p — 1, that paper offers the approximation 

D*(K) w ( ^— ) (13) 



n 

by asymptotic arguments. According to our results in TheoremEl the approximation ( fT3l is indeed a lower 
bound for the distortion rate function and valid for all possible n's. For Q n)P (C) with p fixed and n ^> p, 
a lower bound of an upper bound on the distortion rate function is given in [8] based on an estimation 
of the minimum distance of a code. It is less robust than the result in Theorem El in that it is neither a 
lower bound nor an upper bound, and it only holds for n ^> p (see Fig. El for an empirical comparison). 

Besides characterizing the rate distortion tradeoff, we are also interested in designing a code to minimize 
distortion for a given code size K. Generally speaking, it is computational complicated to design a code to 
minimize distortion directly. In [5] and [12], a suboptimal design criterion, i.e., maximization the minimum 
distance between codeword pairs, is proposed to reduce computational complexity. Refer this suboptimal 
criterion as max-min criterion. According to our volume formula 0, the same criterion can be verified. 
Let the minimum distance of a code C be 5. Note that the metric balls of radius | and centered at Pi G C 
are disjoint. Then the corresponding distortion is upper bounded by 

D (C) < jKfi (B (5/2)) +p(l-Kfx(B (5/2))) . (14) 

Apply the volume formula Q. An elementary calculation shows that the first derivative of the upper 
bound is negative when 



5 < 



A(3p 2 (n — p) 
2 + (3p (n — p) 
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This property implies the upper bound (fLU) is a decreasing function of 5 when 5 is small enough. Thus, 
max-min criterion is an appropriate design criterion to obtain codes with small distortion. Since this 
criterion only requires to calculate the distance between codeword pairs, the computational complexity is 
less than that of designing a code to minimize the distortion directly. 

Fig. 12 compares the simulated distortion rate function (the plus markers) with its lower bound (the 
dashed lines) and upper bound (the solid lines) in (fTTT) . To simulate the distortion rate function, we use 
the max-min criterion to design codes and use the minimum distortion of the designed codes as the 
distortion rate function. Simulation results show that the bounds in (fTTT) hold for large K. When K is 
relatively small, the formula (fTTT) can serve as good approximations to the distortion rate function as well. 
Simulations also verify the previous discussion on the difference between the two bounds. The difference 
between the bounds is small and it becomes smaller as n increases. In addition, we compare our bounds 
with the approximation (the "x" markers) derived in [8]. Simulations show that the approximation in [8] 
is neither an upper bound nor a lower bound. It works for the case that n — 10 and p = 2 but doesn't 
work when n < 8 and p = 2. As a comparison, the bounds (fTTT) derived in this paper hold for arbitrary 
n and p. 



C. Plane Matching: the Rate Distortion Tradeoff 

For completeness, this subsection summarizes the results about the rate distortion tradeoff for plane 
matching. The corresponding proofs follow those for the quantization problem. 

Let a code C £ Q n , P (L) . The plane matching problem is to choose a plane in C to match a random 
plane Q £ Q n>q (L). Without loss of generality, we assume 1 < p < q < n. Denote the size of the code C 
by K. When K is sufficient large, the distortion rate function is bounded by 

(3p (n - q) 



Pp(n-q) + 2 
2T 



c n , P , q ,fsK)-^n (1 + (1)) < D* (K) 



When the required distortion D is sufficiently small, the rate distortion function is bounded by 

pp{n — q) 

1 | (3p(n — q) 



Cn,p,q,P y2r ^ 



0p{n-q) 



D] (l + o(l)) <K*(D) 



, , fiy(.n-q) 

1 ( f3p(n-q) + 2 



Cn, P ,q,f3 V PP{n-q) J 

We finally detail the above (1 + o (1)) errors, and so those in (fTTT) and (fT2ll as well with the following 
theorem. 

Theorem 3: Let a be an arbitrary real number such that < a < 1 and K be sufficiently large 

2 

{{cn, P ,q,i3K) a>(™-9) < 1 necessarily). If L = R and q = p, then 
p (n — p) 



p {n — p) + 2 

2r 



(cn,p, P ,xK)-^ (l - {c ml K)-^~) n - p < D* (K) 



< — P ? 7 (Kc n:P:P:1 ) p("-p) +pexp (- (Kcn^xf a ) 
p(n — p) 
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If L = R and q = p + 1, orL = C and q = p, then 



/3p (n - g) + 2 

2r x 

/3p (n — g) 

If L = R and g > p + 1, orL = C and q > p, then 

(3p (n — q) 



l3p{n-q) + 2 



(c ml3 Kym^n < d* (K) 



2T { -5-7^ r ) , . /3p( 9 -p+l)-2p 



/3p( 9 -p+l) 
2a \ 5 P 



1-0/-, / /= a , \ 2 



(3p (n — q) 

+ pexp (- {Kc n>p>q ^) L ~ a (l - (Kc n>Pt g 

Proof: The proof is given in Appendix |H] and O 
The lower and upper bounds are asymptotically identical. Let p and q be fixed. Let n and the code 
rate log 2 i^ approach to infinity simultaneously with f = lim log2 — . If the normalized code rate f is 

(n,K)— >oo 

sufficiently large, then 

lim D* (K) = PT& . 

(n,K)— *+oo 

On the other hand, if the required distortion D is sufficiently small, then the minimum code size required 
to achieve that distortion satisfies 

lim iog 2 ^P) = ^ 1 (py 

V. An Application to MIMO Systems with Finite Rate Channel State Feedback 

As an application of the derived quantization bounds on the Grassmann manifold, this section discusses 
the information theoretical benefit of finite-rate channel-state feedback for MIMO systems using power 
on/off strategy. In particular, we show that the benefit of the channel state feedback can be accurately 
characterized by the distortion of a quantization on the Grassmann manifold. 

The effect of finite-rate feedback on MIMO systems using power on/off strategy has been widely studied. 
MIMO systems with only one on-beam are discussed in [4] and [5], where the beamforming codebook 
design criterion and performance analysis are derived by geometric arguments in the Grassmann manifold 
Gn,i (C). MIMO systems with multiple on-beams are considered in [8], [13]— [16]. Criteria to select the 
beamforming matrix are developed in [13] and [14]. The signal-to-noise ratio (SNR) loss due to quantized 
beamforming is discussed in [8]. The corresponding analysis is based on Barg's formula dTOb and only valid 
for MIMO systems with asymptotically large number of transmit antennas. The effect of beamforming 
quantization on information rate is investigated in [15] and [16]. The loss in information rate is quantified 
for high SNR region in [15]. That analysis is based on an approximation of the logdet function in the high 
SNR region and a metric on the Grassmann manifold other than the chordal distance. In [16], a formula to 
calculate the information rate for all SNR regimes is proposed by letting the numbers of transmit antennas, 
receive antennas and feedback rate approach infinity simultaneously. But this formula overestimates the 
performance in general. 

The basic model of a wireless communication system with L T transmit antennas, L R receive antennas 
and finite-rate channel state feedback is given in Fig. The information bit stream is encoded into the 
Gaussian signal vector X G C sxl and then multiplied by the beamforming matrix P 6 (£ l txs t0 g enera t e 
the transmitted signal T = PX, where s is the dimension of the signal X satisfying 1 < s < Lt and 
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P on l 8 where the P on 



the beamforming matrix P satisfies ptp = I s . In power on/off strategy, E [XX 1 "] 
constant denotes the on-power. Assume that the channel H is Rayleigh flat fading, i.e., the entries of H 
are independent and identically distributed (i.i.d.) circularly symmetric complex Gaussian variables with 
zero mean and unit variance (CAf (0,1)) and H is i.i.d. for each channel use. Let Y e C Lrx1 be the 
received signal and W e C LrX1 be the Gaussian noise, then 



^L T xs 



Y = HPX + W, 

where E [WW*] = I Lr . We also assume that there is a beamforming codebook B = {P« G 
PjPj = I s | declared to both the transmitter and the receiver before the transmission. At the beginning 
of each channel use, the channel state H is perfectly estimated at the receiver. A message, which 
is a function of the channel state, is sent back to the transmitter through a feedback channel. The 
feedback is error-free and rate limited. According to the channel state feedback, the transmitter chooses 
an appropriate beamforming matrix Pj e B. Let the feedback rate be i?fbbits/channel use. Then the size 
of the beamforming codebook \B\ < 2 Rih . The feedback function is a mapping from the set of channel 
state into the beamforming matrix index set, (p : {H} — > {i : 1 < % < \B\}. This section will quantify 
the corresponding information rate 



I = max max E 

B:\B\<2 R lb <P 



log 



I Lr + P on HP^(H)P; (H )H 



where P on — p/s and p is the average received SNR. 

Before discussing the finite-rate feedback case, we consider the case that the transmitter has full 
knowledge of the channel state H. In this setting, the optimal beamforming matrix is given by P opt = V s 
where V s e C LtXS is the matrix composed by the right singular vectors of H corresponding to the largest 
s singular values [6]. The corresponding information rate is 



2o P t — E H 



^ln(l + P 01 A) 



(15) 



where Aj is the i th largest eigenvalue of HH*. In [6, Section III], we derive an asymptotic formula to 
approximate a quantity of the form E H Ei=i m (1 + c ^i)] where c > is a constant. Apply the asymptotic 
formula in [6]. X op t can be well approximated. 

The effect of finite-rate feedback can be characterized by the quantization bounds in the Grassmann 
manifold. For finite-rate feedback, we define a suboptimal feedback function 



i = if (H) 4 arg min d 2 c (V (P.) , V (V s )) , 

l<i<\B\ 



(16) 



where V (P«) and V (V s ) are the planes in the Gl t ,s (C) generated by Pj and V s respectively. In [6], 
we showed that this feedback function is asymptotically optimal as Rr> — > +oo and near optimal when 
.Rfb < +oo. With this feedback function and assuming that the feedback rate is large, it has also been 
shown in [6] that 



E 



H 



5^1n(l +^supPonAj) 



(17) 



where 



Vsup 



= 1-- inf E Vs 

SS:|B|<2- R fb 

= 1 - -D* (2 Rfb ) . 



min <^(P(Pi),P(V 8 )) 

l<i<|B| 



(18) 
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Thus, the difference between perfect beamforming case < TT3T > and finite-rate feedback case (flTT ) is quan- 
tified by ?7 SU p, which depends on the distortion rate function on the Ql t ,s (C). Substituting quantization 
bounds (ITTb into (fTSt and applying the asymptotic formula in [6] for E H Ki=i m (1 + c ^i)] produce 
approximations of the information rate 1 as a function of the feedback rate Rr>- 

Simulations verify the above approximations. Let m = min (Lt, Lr). Fig. |4] compares the simulated 
information rate (circles) and approximations as functions of R^/m 2 . The information rate approximated 
by the lower bound (solid lines) and the upper bound (dotted lines) in (fTTT) are presented. The simulation 
results show that the performances approximated by the bounds (fTTT) match the actual performance almost 
perfectly. As a comparison, the approximation proposed in [16], [17], which is based on asymptotic 
analysis and Gaussian approximation, overestimates the information rate. Furthermore, we compare the 
simulated information rate and the approximations for a large range of SNRs in Fig. |5] Without loss of 
generality, we only present the lower bound in (fTTT) because it corresponds to the random codes and can 
be achieved by appropriate code design. Fig. 0a) shows that the difference between the simulated and 
approximated information rate is almost unnoticeable. To make the performance difference clearer, Fig. 
|5fb) gives the relative performance as the ratio of the considered performance and the capacity of a 4 x 2 
MIMO achieved by water filling power control. The difference in relative performance is also small for 
all SNR regimes. 



VI. Conclusion 

This paper considers the quantization problem on the Grassmann manifold. Based on an explicit volume 
formula for a metric ball in the Q njP (L), sphere packing bounds are obtained and the rate distortion 
tradeoff is accurately characterized by establishing bounds on the distortion function. Simulations verify 
the developed results. As an application of the derived quantization bounds, the information rate of a 
MIMO system with finite-rate channel-state feedback and power on/off strategy is accurately quantified 
for the first time. 



Appendix 

A. Proof of Theorem [7] 

The proof is divided into three parts, in which we calculate the volume formula for the 1 < p < q < ^, 

§ < P < q < n and l<p<^<q<n cases respectively. 

1) 1 < .P < <? < § case: First we prove the basic form 

H (B (5)) = c^rff**-* (l + + o (5 2 ) 

for 5 < 1. Afterward, we calculate the constants c„ )3 , )?)/ g and c^ pqf3 . 
The volume of a metric ball /i (B (5)) is given by 

MS (<*))= due, (19) 



where the differential form dfie is the joint density of #j's. For convenience, we introduce the following 
notations. Define r,i = cos 6^ and order n's such that r { < rj (9{ > 8j) if i < j. Define r = [r 1; • ■ • , r p ] 
and also 



i A p( r2 )i=n 



2 _ 21 

r- r i | 



i<j 
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Recall that (3 = 1 for L = R and (3 = 2 for L = C. With these notations, the invariant measure d^e can 
be written as follows [11]. 

dp, = v mP |Ap (r 2 ) f f[ {(rf)^- p+1) - X (l - r 2 ) f ("-p^+iM drf ) > (20 ) 



where the constant v n;P;qi p is given by 

v A r(i + f)r(j( B -« + i)) 

<w 11 r (f i + i) r (f („ - P - i + 1)) r (f ( ? - > + 1)) ■ (/ " 

To get the form Q, we perform the variable change S 2 Xi = 1 — rf for 1 < % < p. Under this 
transformation, the integral domain 

Ds ~ { r : i 1 - r i) - ^' < r 2 < 1 for z = 1, • ■ • ,p| 

is changed to 

D 1 = jx : y^Xj < 1, < Xj < for % = 1, • - ,p| 

= |x : Xj < 1, < Xj for % = 1, • ■ ■ ,p| , 
where the last equation holds since 5 < 1. Thus 

(5 (5)) = /••■/|A(x)fn((l- dx, 

Next note that 

(1 - S'x^-^- 1 = 1 - U (q - p + 1) - 1 J x,5 2 + o (5 2 ) , 
and so we are able to express the volume of B (5) in the desired form with 

w y • • • y ia (x)i' 3 n ((x,)^-^-^ 1 )- 1 <&< 



and 



/•-•MiA(x)i^nr=i ((^) |(n - p - ?+1) - 1 ^ 



2 / /•••/iA(x)i /3 m = ife) l(n ^ +1M ^ 



In order to calculate the constants Cn,p, q ,p and c^ pqf3 , we need the following lemma [18]. 
Lemma 1: It holds that 

7-1 



r(7) 



i=i / i=i 

m 



T (7 + m + + |p (p — 

A r(a + |(p-0)r(i + fQ 

x ll f(i+I) 
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where < m < p, 3ft (a) > 0, 3ft (7) > 0, 3ft (f ) > - min Q|, and the integral is taken over 

Proof: This is Selberg's second generalization of the beta integral. See [18, Section 17.10] for a 
detailed proof. ■ 
According to Lemma [H we have that 

j... j iA( X )rn((x i )^"- p - 9+i) - 1 ^) 



i=l 



1 



r (p§ (n - p - g + 1) + §p (p - 1) + 1) 

' r(l + |*)r(f (n-p-g + l) + f (p-0) 



x 



n 



r(i + f) 

n r(l + f.)r(f (n-q-j + l)) 

r(fp(n-g) + l)fj r(i + f) 

Substituting the formula (I2TT) for f n ,p,g,/3 into the integral expansion of c n ^ p ^ yields 

1 *L r(f(n-i + l)) 



Cn,p,q,f3 



r(fp(n-g)+l)n r (| (g _ i + i)) 



after some simplifications. 

The constant ^ can be calculated in a similar way. Lemma [T] implies that 

r~t i=l 



f (n - p - g + 1) + f (p - 1) X T (1 + § i) T (I (n - g - % + 1)) 



n 



Therefore, 



r(fp(n-g) + 2) 1A T (l + f ) 

f("-g) f 1 A r(i + fQr(f(n-. 
| p ( n _ g ) + i ^r(| p ( n _g) + i) 11 r(i + f) 



which completes the proof of the 1 < p < q < | case. 

2) I < p < g < n case: This computation is closely related to that for the 1 < p < g < | case. 

To see the connection between the | < p < g < n and 1 < p < g < ~ cases, we define the generator 
matrix and the orthogonal complement plane. For any given plane P E ^ n p (L), the generator matrix 
P g L nxp is the matrix whose p columns are orthonormal and expand the plane P. The generator matrix 
is not uniquely defined. However, the chordal distance between P E Q np (L) and Q E G n ,q (L) can be 
uniquely defined by their generator matrices. Indeed, 

d\ (P, Q) = min (p, g) - tr (P+QQ+P) , 

where P and Q are generator matrices for the plane P and Q respectively. It can be shown that the 
chordal distance is independent of the choice of the generator matrices. The orthogonal complement 
plane is defined as follows. For any given plane P E Q rhP (L), its orthogonal complement plane P L is the 
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plane in G n ,n- P (L) such that the minimum principle angle between P and P x is |. It is straightforward 
that P^P- 1 = where P and P^ are the generator matrices for P and P 1 - respectively, and the matrix 
is the p x (n — p) matrix with all elements 0. 

With the definition of the orthogonal complement plane, the chordal distance between P and Q can be 
related to that between P L and Q ± . The relationship is given in the following lemma. 

Lemma 2: For any given planes P e Q n>p (L) and Q e Q n>q {h), let P 1 - e Q n ,n-p (L) and Q -1 G 
Q n>n -q (L) be their orthogonal complement planes respectively. Then 

d 2 c (P,Q) = d 2 c (P ± ,Q ± ). 
Proof: This lemma can be proved by the generator matrices. Let P, Q, P^ and Q- 1 be the generator 
matrices for P, Q, P L and Q 1 - respectively. Without loss of generality, we also assume that 1 < p < q < n. 
Then 

p = tr(pt[Q|Q^ [QIQ^P) 

= tr (P^Q+P) + tr (ptQ ± (Q- 1 ) 1 p) , 
where the matrix [Q | Q- 1 ] is the one composed of Q and Q- 1 . Similarly, 
n-q = tr^Q^plP- 1 ] [P^V) 

= tr ^(Q ± ) t PP t Q ± ) + tr ((Q^ P x (P^) 1 Q x ) . 

Then 

d 2 c (P,Q) ( = } p-tr(ptQQtp) 

= tr(ptQ^(Q^) f p) 

= n-g-tr((Q^) t P^(P^) t Q^) 

where (a) and (b) are from the definition of the chordal distance and the facts that min (p, q) = p and 
min (n — p,n — q) = n — q. ■ 
By this lemma, the connection between the | < p < q < n case and 1 < p < q < | case is clear. The 
volume formula for the | < p < q < n case can be calculated as follows. 

H{B P {8)) = Pr {Q e g n , q (h) : d 2 c (P,Q)<8 2 ) 

= Pr (g X G S n ,n- 9 (L) : d\ (P x , Q 1 ) < 5 2 ) 

where Bp (5) and B P ± (5) are the metric balls in Q nA (L) and Q n , n -q (L) respectively. Therefore, the 
results for the 1 < p < q < | case can be directly applied by letting p' = n — q and q' = n — p. Finally 
after some simplification, we have 

H (Bp (6)) = c mP 5^ (l + c^5 2 + o (S 2 )) , 

where 

1 yf r(f(n-» + l)) 

Cn ' P ' q ' P r(lp(n- q ) + i)l\r(l(n-p-i + i)Y 

and 

(i) (P, , n ^ IP""? 
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3) l<p<T;<q<n case: This computation is again related to that for the 1 < p < q < | case. 

Similar to the | < p < q < n case, the connection between the | < p < q < n case and 1 < p < q < | 
case can be revealed by the generator matrix and the orthogonal complement plane. Let p' = min (p, n — q) 
and q' = max (p, n — q). Then 1 < p' < q' < |. For any given planes P 6 p (L) and Q E Q nA (L), let 
Q 1 " £ {? n , n - g (L) be the orthogonal complement plane of Q. Let P, Q and be the generator matrices 
for P, Q and Q x . Then 

d 2 c (P,Q) = p-tr(ptQQtp) 

= tr^Q^Q^p) 
= p'-dl{P } Q^). 

Therefore, 

fi(B P (5)) = Pr(Qeg n , q (h): d 2 c (P,Q)<5 2 ) 

= Pr (Q X G S n , n _ ? (L) : d\ (P, Q L )>p- 5 2 ) . 

Now calculate the volume formula. Note that 

fM (B P (5)) = Pr {Q 1 - E G m (L) : d\ (P, Q x ) > p' - 5 2 ) . 

Then 

V (Bp (5)) = J "J d ^ e <p'-i' 

o<^ 2 <i,E?1i(i-'- s 2 )>p'-<5 2 

d^e, P ',q'i 

where dfi ep i q i is the invariant measure with parameter n, p' and q'. Substitute the form for dfie^/^i d^Ut 
into the above formula. Then 



//(Bp (*)) = t;, 



c^,^ 9 ' (l + ? , / + o (5 2 ) 
where the first equation comes from the variable changes 5 2 Xi = r 2 (1 < % < p'), v np i q i p is defined in 

<ED, 

, r 

D\ = < x : < 1, < Xi for i = 1, • • • ,p' 



C n ,p>,q',P = V n ,p',q',l3 J "J I A ( X ) 1^ fl ^ P + ^ ' 



and 



P ',q',p — y o v° v y ~ / / /~r 
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Applying Lemma [T] and after some simplification, we have that 



1 fr r(f(n-i + l)) 

Cny '^ r(f P Y + i)Mr(f(n-p'-i + i)) 



and 

c n.»'. './3 = - 77 ( n ~ v - q + i) - i 1 - 



Summarily, if 1 < p < | < q < n, 

» (B (5)) = c n , m ^-* (l + c« ?) / + o (5 2 ) 



where 



and 



l -prn-g r(f(n-i+l)) 

C„,p, g ,/J <j ^ r(|(n-i+l)) 

r(fj.(n-,)+i) 1 A<=i r(§ (n-p-i+i)) 11 p + q - n 



5. Proof of the lower bound on D* (K) 

Assume a source Q is uniformly distributed in £/ np (IL). For any codebook C, define the empirical 
cumulative distribution function as 

i^ jC (x)=Pr {Q: fmm^(P,Q)J < *} . 

Then the distortion associated with the codebook C is given by 

D(C) = J x-dF^dx). (22) 

The following theorem gives the empirical distribution to minimize the distortion. 
Lemma 3: The empirical distribution function minimizing the distortion for a given K is 

if x < 
F d % fi (x)=\ K-n{B(y/x)) if , 

1 if x* < x 



where x* satisfies K ■ fi (B (y/x*)) = 1. 

Proof: For any empirical distribution c (x), 



F dlc (x) = Pr<{Q: ( mm d 2 c (P, Q) ) < x 



= Pt(uK =1 {Q: d 2 c (P h Q)<x}) 

K 

< £>r{Q: d 2 c (P i ,Q)<x} 
i=i 

= K.fA(B(y/Z)). 



Thus 

F dlc (x)<mm(l,K-fi(B( v ^))). (23) 
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Therefore, 



o 



(«) 



p 

x ■ dF d 2 C (x) — I x ■ dFj2 C (x) 

r-p 



rp rp 

/ F dlc( x ) dx ~ / F d? c ,c{x)dx 
Jo Jo 

/ (K ■ fi (B (v^)) - F Alfi (x)) dx+ (1 - F dlc (x)) dx 

Jo Jxo 



(b) 

>o, 

where (a) follows from integration by parts, and (b) follows from d2"31) . 
From Lemma it is clear that 

CP 



D*(K)> f x-dF* lfi {x) 
Jo 
rp 

= / x-dKF{x), (24) 
Jo 



where F (x) = fi (B (yfx)). 



1) Proof of Theorem^}; Theorem|2]is proved by substituting the volume formula © into (|24h . Another 
way to prove it is to apply the lower bound in Corollary El whose proof is more involved and given in 
the following. 

2) Proof of the lower bounds in Corollary The difficulty to calculate d24l is that we don't know 
the exact F (x) for some cases. To overcome this difficulty, we construct a further lower bound on (l24l) . 

For all cases except the /3 = 1 and q = p case, a lower bound on d2"4l is constructed as follows. Let 

F (x) = c njm ,pxf p( - n -d and x satisfy KF (x ) = 1. Since F (x) < F (x) (Corollary H), KF (x ) < 
KF (x ) = 1. But KF (x*) = 1. We have x < x*. Therefore, 

X* 

x ■ dKF (x) 



> 



/•X* 

/ {l-KF{x))dx 
Jo 

rxo 

/ (l-KF{x))dx 
Jo 



-co 



> / (1- KF (x))dx 



xo 



x ■ dKF (x 
Pp (n — q) 



(c n , P , q ,pK) m^i. 



(3p(n-q)+2 

For the case (3 = 1 and q = p, the computation is more complicated. The following lemma is helpful. 

Lemma 4: Let a = \p{n — p). Let F (x) = c n:P:P: iX a and x satisfy KF (x ) = 1. Let F ub = 
c n , P ,p,ix a (1 - and x ub satisfy KF ub (x ub ) = 1. Let F ubub = c n ^ 1 x a (1 - xo)^ and x ubub satisfy 
KF ubub (x ubub ) = 1. Then 

/ x ■ dKF ubub (x)< x- dKF ub (x) < x- dKF (x) . 
Jo Jo Jo 
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Proof: Similar to the arguments for all the cases except the (3 — 1 and q = p case, it can be proved that 

x u b < x* < x and f* ub x ■ dKF ub (x) < x ■ dKF (x). Then (1 - x)~^ < (1 - x )~^ for x G [0, x ub \. 
It implies F ub (x) < F ubub (x) for x E [0,x ub \. Therefore, x ubub < x ub and f* ubub x ■ dKF ubub (x) < 
J* ub x-dKF ub {x). m 

We calculate J* ubub x ■ dKF ubub (x) as follows. x ubub = x (1 — x )~ . 

r^ubub 

I x ■ dKF ubub (x) 
Jo 

^^ubub 



= / x-dKF (x) (1 -x yz 
Jo 



p(n 


-P) 


p(n- 


p) + 2 


p (n 


-p) 


p{n — 


p) + 2 


p (n 


-p) 


p{n — 


p) + 2 



Kc n>PtP>1 (l- x ) 2 x 



\p(n-p)+l 



ubub 



Xubub 
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C. Proof of the upper bound on D* (K) 

To get an upper bound on D* (K), we shall compute the average distortion of the random codes. Let 
Cand = {-Pi, ' ' ' , Pk) be a random code whose codewords Pj's are independently drawn from the uniform 
distribution on Q n v (L). For any given element Q € Q n:P (L), define X^ = d 2 c (Pi, Q), 1 < i < K. Then 
Xi's are independent and identically distributed (i.i.d.) random variables with distribution function 

F (x) = // (B (v^)) = c m px L * (1 + (a;)) . 

Define Wk = min (Xi, ■ ■ ■ , Xk). Then 

E C ran d [ D (Crand)] 



E C] , 



min dUQ,^) 
min t^(g,Pi) 



P% ^C ra nd 

= E Q [E^ [W k ]] . 

To calculate Ey^ [Wk], we need to know the distribution of Wk- To derive it, the the following lemma 
is useful. 

Lemma 5: Let JQ's 1 < i < K be i.i.d. random variables with distribution function F (x). Let W K = 
min (Xi, ■ ■ ■ , X^). Then 

exp(-KF(a;)) > Pr (W K > x) — (1 — F (x)) K 

where the upper bound holds for all x. 

Proof: See [19, page 10]. ■ 
With the above upper bound on the distribution function of Wk, we derive an upper bound on E Wk [Wk]- 
In the following, we use E [•] instead of E Wk [•] for simplicity. Let F ib (x) be an arbitrary distribution 
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function such that F ib (x) < F (x). It is clear that Fi b (x) is zero if x < 0. Then 

rp 

E [W K ] = / Pr (W K > x) dx 
Jo 

(a) fP 

< / exp(-KF (x))dx 

Jo 

rp 

< / exp (—KFib (x)) dx 
Jo 

px rp 

< / exp (— KFi b (x)) dx + / exp (— KF\ h (x )) dx 

Jo Jx 



< / exp (—KFib 0)) dx + pexp {-KF lb (x )) , 
Jo 

2 

where (a) follows Lemma |5J Here and throughout, xq = (c n ^ tq ^K) ~ a , a G (0,1) is an arbitrary 
real number and K is large enough to guarantee x < 1. 

If 1) L = R and p < q < p + I, or2)L = C and q = p, then we can take F w (x) = c n:P:(]t pX2 p ( n ~ qS> 
(Corollary [l]and 0). We have 



E [W K ] < J exp (-Kcn^px dx 

+ pexp ( -Kc n)P) q )/3 ((Kc n)P) q )/3 )~~ W"- 9 ) a j 



0p(n-q) 



= (Kcn^p) M^T) I e y ■ dy* 



1 



+ pexp (- (Kcn^p) 1 a ) 

2 



2T 



2 

where (b) is from the variable change y = Kc n ^^x^ n -^ . Thus, for any given n, p and q, 

2T 



lim E 

K->oo 



K (3p{n-q) 



and so, when K is large enough, 

2T 



D* (K) < - ^ n ~ q) ' ( Cn ^Kyw^ (l + o (1)) . 



f3p (n — q) 

If 1) L = M and q > p + 1, or 2) L = C and q > p, then F (re) > c nvq {ix a (1 — x) 7 where 



a = | p ( n - q ) and 7 = Mi-rt-i) - p (Corollary [U and Q. Note that (1 - i) 7 > (1 - ,t„) 7 for all 
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x E [0, x ]. We take F lb (x) = c niPiq:f3 x a (1 - x o y . Then 

E [W K ] < / exp (-Kc ml3 x a (1 - x ) 7 ) dx 
Jo 

+ pexp (-irc niPi g i/3 Xo (1 - x o y) dx 

( £ ) _i / f°° i 

< (Kcn^p) <* [1 - (Kcnj,^) a J J e v -dy* 

pexp {- {Kc n ^f~ a (l - (Kc n ^y«y^j 

= ^ (*<w)~* (l - (Kc m ^Y" 

+ pexp (- {Kc n ^ p ) 1 ~ a (l - {Kc n ^)~^y^j , 
where (c) follows the variable change y = Kc njPjq> px a (1 — x o y. Once more, for any given n, p and q, 

2T (/8p(n-5)) 



lim E 



pp (n - q) 



(Cn,p,q,f3) , 



and 



2r 



L>* (K) < 



/3p (n - g) 



^p,^) (1 + 0(1)), 



for sufficiently large if. 
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Real Grassmann Manifold 




(a) Real Grassmann manifolds 
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(b) Complex Grassmann manifolds 
Fig. 1. The volume of a metric ball in the Grassmann manifold 



distortion rate function for G (C) 




Fig. 3. System model for a MIMO system 
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Fig. 4. Performance of a constant number of on-beams v.s. feedback rate 
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5. Performance of a constant number of on-beams v.s. SNR 



